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ABSTRACT 


Pi xcdweomrryeiuravion feedback control oftenmetomissuic 
Pio nom Cie Gesifl of A linear control systemaeince 
design freedom is limited to the adjustment of the free 
parameters of the system. Analytical methods prove to be 
valuable in these circumstances and the aim is the optimi- 
Zation of a selected index ei yperformance. A Saeeestea 
index is developed for a regulator system from a comparitse 
of the dynamics of the subject sy&Stem witha those of a 
desired model along a transient state space trajectory. 
The result is an Optimal system dependent on the into 
state chosen for the transient response. The effects of 
the initial state on the optimal system are investigated 
and tne final selection is based on statistical consider 


a 1 Oms . 
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CHAPTER I 


INTRODUCTION 


Tne design of feedback control systems sometimes in- 
volves the selection of free parameters in the plant and/or 
compensator in order to satisfy the requirements of the per- 
formance specifications. The configuration of the system is 
assuned to be fixed, i1.e., a compensator has already been 
chosen and the design freedom is limited to the adjustment 
of cne ‘or more parameters, fhe trial amd error precess or 
the selection of several parameters represents a very labor- 
jous and sometimes impossible task for the classical design 
by root locus methods. Analytical methods are appropriate 
Mmmenia case, <nehsureor QUMGGe Jewestanmieeee Moaslané 
system in question by the construction of an index of per- 
formance, which upon minimization (or maximization) will 
yield tne optimum values for the free parameters of the 
systen. 

Different criteria could be adopted for the construc- 
tion of tne index of performance, all normally based on 
physical reasons and each satisfying particular requirements. 
Many performance indices originating from good engineering 
reasoning lead to analytically insolvable problems. A com- 
promise in the selection of the index is then necessary to 
obtain a solution. The classical forms of the performance 
specifications, such as rise time, maximum overshoot, 
settling time, etc., in the tite domain or bandwidth, phase 
and gain margins, etc., in the frequency domain, are trans- 


formed into a single measure by conceiving an ideal model 


( 


whose behavior meets all the performance specifications. A 
mathematical relationship established between some variables 
of this model and the corresponding variables in the actual 
system serves as a measure of quality as the free parameters 
are changed in a search for the optimum adjustment. The 
time integral of the difference-Ssquared between the output 
of the model and the output of an actual system is a practi- 
cal index of performance when the system is subjected to 
transient input signals. Here the model has a transfer 
function equal to unity. This is usually solved with the 
aid of Parseval's re ose ee tle 

An index of performance has been derived for the regu- 
lator ero This index of performance, or cost func~ 


tion, is a quadratic function of the state vector, 


J - fx"Qxut (1) 


The matrix of the quadratic form, Q, is the weighting 
matrix and contains the free and fixed parameters of the 
actual system and some weighting factors determined by the 
model or ideal system. The model obtained from the per- 
formance specifications is of the same order as the actual 
system. The weighting factors are adjusted in such a way 
that if all the parameters of the actual system were free, 
the cost function would exhibit an absolute minimum at the 
model values. If the actual system is constrained to have 
less than n independent free parameters, a minimization of 


the cost function will Yield the optimum values for these 


parameters, 

The measure of quality is obtained by disturbing the 
system with some initial conditions. For a stable systen, 
the state space trajectory ends up at the origin for large 
values of time. This type of excitation is characteristic 
Of. wbhewT een batomsDGOD keel nimizine this, cost, fonction 
the trajectory of the subject system is forced to approximate 
as closely as possible the trajectory of the model. If all 
n parameters of the system were free, it would indeed simu- 
late identically the dynamics of the model and hence satisfy 
in full agreement the performance specifications. In any 
etner case, the model behavior will not be reproduced, but 
the actual system will attain its best performance as 
measured by this optimization ¢riterion. 

The importance of this method over the common types 
of indices is the consideration given to every part of the 
dynamics of the system and not to a single manifestation as, 
hear example, the, error in position... The contributions from 
all the state variables are properly weighted according to 
the characteristics of the model, which carries the infor- 
mation of the performance specifications, 

The evaluation of the integral in equation (1) is done 
by Liapunov functions with vector-matrix methods, and the 
final form is shown to be a quadratic function of the 


initial state, 


J = x1(0)Px(0) (2) 


The matrix of this quadratic form, P, is related to the 


Weighting matrix and the actual system's matrix. 

The nature of this solution implies that the cost 
function will be different for different initial states and 
hence the optimum system will depend on the initial condi- 
tions chosen for the calculation, 

In this paper, a similar index of performance is de- 
veloped, based upon a different approach than that described 
above. To obtain a measure of performance, at every point 
on the state space trajectory of the actual system, a com- 
parison is made between the dynamics of this system and the 
dynamics the model would have if it were at that particular 
porns 

The index of performance can be expressed in the same 
form as equation (1) and:‘hence is also a function of the 
initial state. However, the weighting matrix is now direct- 
ly obtained from the coefficients of the model and actual 
system characteristic equations. This is a much more direct 
method than that for obtaining the Q matrix in equation (1), 
which relies on complicated algorithms for the weighting 
factors. In the analysis of this index it is found that 
the results are exactly the same as for the index of equation 
(1). Further investigation shows that both cost functions 
are intimately related by a constant factor and an additive 
constant. The effects on the optimum system of the initial 
state chosen for the calculation are investigated in this 
paper and the final decision on the design is based on the 
statistical characteristics of the disturbing initial con- 


ditions. 
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The new concept developed also suggests a time-varying 
optimization scheme by which the actual system would be made 
to follow exactly the state space trajectory of the model by 
continuously varying its free parameters so as to match at 
every point the dynamics of both systems. Examples are pre- 


sented to complement the theory. 


Al 


CHAPTER IT 


THE INDEX OF PERFORMANCE 


The first step in the process of analytical design is 
to establish the concept on which the index of performance 
Will be based. The objective of this performance index is 
to include in a single number a measure of quality for the 
performance of the system. Factors such as mathematical 
difficulties will perhaps force the designer to modify his 
performance index in favor of mathematical simplicity, hope- 
fully keeping the essence of the idea that originally moti- 
vated it. Once the configuration of the system is fixed and 
the index of performance has been defined, the design 
proceeds in a straightforward manner. For the type of index 
to be considered here, the performance specifications are 
transformed into an ideal or model system thus fixing the 
location on the s-plane of the required zeros and poles of 
the transfer function. This model is of the same order as 
the actual system and fully satisfies the specifications. 
The model and actual systems can be represented by the signal 
flow graphs shown in Figure 1. 

The feedback coefficients of the model are defined by 
the performance specifications. Some of the coefficients 
of the actual system are constrained by the physical nature 
of the plant, and the others are functions of the free para- 
meters. If all coefficients in the actual system were free 
it would be trivial to reproduce the model. Generally this 


is not the case and the free parameters must be adjusted 


ie 


according to seme ciitveri one 
A measure of quality which originates from a good 


physical concept is given by 





J = [[t- ae] # o 


where the integrand is the sguare of the difference between 
the nth derivatives of the model and actual systems differ- 
ential equations respectively, when the system is excited by 
initial conditions only. This index gives a comparison of 
the complete dynamics of both systems in terms of the high- 
est derivatives, and does not require any weighting factors 
Since it implies a direct relationship between tne two 
systems, The minimization of equation (3) as function of 
the free parameters results in an actual system whose dynamic 
behavior is a close approximation to that of the model. 
However, tne development of this performance index into a 
practical form presents the impediment of requiring the 
solution of both system differential equations, which shadows 
the merit of its physical meaning and renders the integra- 
tion of equation (3) very difficult. [It is then necessary 
to conceive a similar structure for the index of performance 
but one which permits convenient mathematical manipulations, 
A Similar way of evaluating the performance of tne 
system can be obtained from a comparison along a trajectory, 
of the nth derivative of the actual system differential 


equation with that of the model when both systems are in the 
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same state. It is again assumed that the system is driven 
by initial conditions only. Tne motivation for this index 
arises from a consideration of a second order process. Con- 
sider the system shown in Figure 2. The system differential 


equation is 
rane AX + a,X = 0 (4) 


where a4 is the free parameter 
The performance specifications require a model system char- 


acterized by the differential equation 
¥ + BY + &yY = 0 (5) 


The actual system trajectory in the phase plane is shown in 
Figure 3, At any point on this path, x is determined by the 
coordinates of the phase plane and the coefficients of the 


differential equation 


x, = ~2nX, -a4X4 (6) 


Also, equation (5) determines the value of y at any point on 


the phase plane 
Y, = “BV, -8194 (7) 


where the subscript 1 indicates the particular value of the 
variables at any given point. 

Along any trajectory of the actual system there will be 
a continuous function of the state giving the value of y. 


THis Function. te 
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model actual system 


Fig. 1. Signal flow graphs for model and actual system 





Fig. 2. A typical second order system 
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The ontimization criterion is the minimization of the differ- 
ence squared between the variables X and y along the actual 


system trajectory. The index of performance becomes 


a = [18 -5re (9) 


when actual system and model are in the same state, Tnen 


' pa ; a 
T=} [lg-%)¥ + (qa) | dt a 


Figure 4 shows a Signal flow graph for the calculation of 
the index. 

In eearenine for the optimum adjustment of Buys tne 
trajectory along which this comparison is made will change 
until the value of the index is a minimum. This will be the 
closest trajectory to that of the model and the performance 
of the actual system will have been optimized according to 
this criterion, In terms of the dynamics of motion, we 
Wish to make the acceleration of the actual system at any 


point on the phase plane a best approximation to the accel- 


eration tne ideal system would have at that point. 


This idea is extended to an neh order system 
CL (ce ae +22 4 ax = 0 (11) 
ltt te at 5 
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[x(0) »x(0)| 





Fig. 3. Phase plane trajectory for system of Fig. 2 





Fig. 4. Signal flow graph showing the structure of 
the index 
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where the index of performance becomes 


” ol 
J -{(¢ [ae real ‘ pee? 


when actual system and model are in the same state, 











Tien 
_ _q.\d"s ONE , 
i 7 (= 20) Spe + fou Soe) Se + 
see 2, (13) 
2 (92. 2) “a (94-44) & | olt 
By defining the state vector 
x ate 
z, is 
2 dt 
x = | - (14) 
{htc 
an au 
this index can be written in quadratic form 
2 one 
us -/ |( Ju-On) 2 a v iat Co x i 
oe 
a (9, = 94) xf 5 8 2 (9n-&n)(9n.s ~2n-1) XX) 
Me ee caaae a) 2% | dt 
(15) 
[cy Xg --- Xy] by Dig >> by] fy 
= o> > aie (16) 
° byt ~ = =e b,.., x 


where byj = (8; - a,) (8; - ae i= 1,2,.....n0and) j =a 
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Then 
4 =| ® B x dt (17) 


where the matrix of the quadratic form, B, is a positive 
semi-definite, symmetric matrix. 

This is the general form of the index of performance 
as given in Equation (1). It is solved by vector-matrix 


methods using Liapunov's Theoreml? J, The solution is 
iE 
J = x (0) P x{(0) (18) 


wnere P is the matrix of the Liapunov function and is re- 
lated to the weighting matrix B and the actual system matrix, 


A by 


2? 


if 
-B = A, leacle PA, Wisi) 


The B matrix is directly obtained from the model and 
actual system coefficients as shown in equation (16). The 
actual system matrix is defined by the matrix differential 


equation in canonical phase variable form 


x= A, x (20) 
wnere 
r 
| Oa. | 
er ce | 
| z 
i, ee | | ~~ Pn ( 1) 
| h ‘ 4 
“QA, -Q, aa liar -Qn 


ouega. tlow eraphns are helpfiub in solving equation 
(bo) for the elements of the Pp matrixL! J, Thewcesil lt Saorue. 


second order system are stated here for illustration, 





a arp ccc 205 W5 
His - by + Tite aie a 

oe | b _ P (220) 
(a 20," 21 

ey he b 
Pez 20> 22 20402 NI 


The cost function, equation (18) becomes 


J = se, (OVE + 2x, (0)X,(0)py2 + X5°(0)Po (23) 


The minimization of the cost function will yield the 
optimum adjustment for the free parameters. Because of the 
nature of this artifice, the optimum system will depend on 
the initial state chosen for the calculation. This problem 


is discussed in Chapter III. 
A TIME-VARYING OPTIMIZATION METHOD 


The minimization of the index of performance developed 
in the last section produces a set of parameters for which 
tne dynamics of the system is a best approximation to that 
of the model. In terms of the state space, this means that 
the actual system trajectory is now the closest trajectory 
to that of the model which can be achieved with this con- 
strained system under the effect of a given set of initial 
conditions. The index is now a measure of the total error 
incurred by the system in going from this particular initial 


state to the origin of the state space. Because it isa 


ZO 


quadratic function, it gives more importance to large errors 
than to small errors but is not capable of reducing the mag- 
nitude of any panticular emnuor imeperformance. 

The criterion on which this index is based suggests a 
possible way to obtain with the actual system an exact re- 
production of the model dynamics if the free parameters are 
allowed to vary continuously according to some function of 
the present state, Actually, it would be enough to have 


th order system to implement 


only one free parameter in ann 
this scheme. 

Consider again the second order system of Figure e2. 
The actual system and model differential equations are given 


in equations (4) and (5) respectively. When both systems are 


in the same state. 


and the condition for the two state space trajectories to 


be the same (and so the dynamics of both systems) is 
j= 


or 


ie 


(oo eg Me Box = a,X + ay 


then 


(24) 


Pa] Ps © 


a =; oy “+ (S575) 


“al 


Bquation (24) is the required functional relation be- 
tween the free parameter and the state of the system. This 
1s plotted in ’Ffigure™S> for aytypical case. 

The curve shows a discontinuity whenever the trajec- 
tory crosses the x-axis in the phase plane. Since infinite 
gain is not practically feasable, the function can be approx- 
imated as shown in Figure 5. A cnange in the sign of the 
gain factor from time to time is also required. The reason 
for the discontinuities is that when x = Q, the free paran- 
ewer: ai? whicn modifies the variable x in the feedback path, 
has no effect on the behavior of the system, now solely con- 
trolled by the variable X and its fixed feedback parameter 
os Thus, an infinite gain is the result of the impossi- 
bility of altering the system's trajecevory. 


th order system, the come 


For the general case of ann 
dition under which the actual system will duplicate the 


model dynamics is 


ned 
9 se OS sic ici _ a 


dt does 
n-1 
U,< ae Qs a i Cy a (25) 


If we assume that a, is the only one free parameter, then 


iL 


d 
a, = a seals = Ce) fat +4 09 285) 


Equation (26) shows that only one free parameter is 





dx /d tr (26) 
oo 


needed to reproducé the model dynaitics if thre state coreume 
system is known. Normally the fact that the state vector 


1s completely measurable implies that the feedback paths 
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can be modified to copy the model characteristics exactly. 
iaisearscussilon Wouldwappinein these special cages in @Phich 
the state vector is measurable but cannot be used directly 


in the feedback channel, 
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.CHAPTER III 


THE OPTIMUM SYSTEM AS A FUNCTION OF THE INITIAL STATE 


The index of performance that will yield the optimum 
Ser is a function of the initial state, We wish to in- 
vestigate the effect of the initial conditions on the 
optimum pole locations. These en titen poles will occupy 
certain regions on the s-plane according to the constraints 
pam tie Sh and one requirements for stability. Also of 
Se * the knowledge of how the extension and shape of 
amen eae ‘vary as a function of the relative position 
of the model poles with respect to the actual system's 

attainable roat regions. 

| | Bane: Mie normalized initial state vector wj1ll be con- 
sidered in this discussion since all other points in the 
state space are proportional to this unit vector and differ 
in cost function only by a constant factor. For each 
orientation of enc art initial state vector there is a 
PanGe pondina optimum set of Roles for the actual system 
and a value of the cost function which is a minimum for 
this particular choice of initial conditions, Any other 
Bape neat ken Will yield a different set of optimum poes. 
a it will be sean thats a of the symmetry 
ee by the index of performance, only a nalf of tne 
state space need ie; 1S searched for optimization. Tne 
design must provide a system whose performance is optimum 
in some modified sense for all possible initial conditions. 


The selection of such a set of optimum parameters must 


25 


be done on the basis of the statistical properties of the 
disturbance. If the probability density function can be 
established, then the average value of the cost incurred by 
the system can be found for all sets of optimum parameters, 
The set that provides the minimum average cost is the 
design selection. This, of course, implies the consider- 
ation of the behavior of the control system as a regulator. 
It is not possible to obtain an analytical relationsnip for 
the optimum poles as a function of tne initial state and 
the proximity of the desired model poles for a given con- 
strained actual system. However, the general trend is that 
the regions containing the optimum poles reduce, and finally 
collapse to a point, when the model poles get closer to and 
reach the permissible regions for the actual system. Simi- 
larly, the same effect is obtained by geradually releasing 
the constraints of the actual system but keeping the same 
characteristics of the desired model. These effects will 
be discussed for the second and third order systems by 


means of examples. 
SECOND ORDER SYSTEM 


The results of the optimization scheme developed in 
Chapter II will now be discussed for a typical second order 
system, Two parameters define the dynamics of this system. 
The method is applicable when one of this parameters, 
usually in the plant, is fixed. With only one degree of 
freedom, the possible poles for this system must lie ona 


definite locus. A pair of required poles outside this locus 


26 


is simply unattainable and the optimization process must 
be-employed . 
Consider the second order system shown in Figure 6, 


‘The transfer function is 


(27) 


IO 


, k 
- ae 
_S +s4+ k 


The free parameter is the gain of tne system wnereas tne 
plant contains the fixed parameter. The performance 
specifications require a model defined by the transfer 


-function 


Ta 
weN 


a 4 (28) 
a eet 7 fo ena 
gra oes ae 


> 
es Pony 
sa Te) 
re iy 


It is noted that the model poles are irate ea dines meee 
the actual system, It is Dian mene to find an adjust- 
Pent for k that will give the best performance with respect 
to the mts according to tne criterion previously establish- 
a - . ; 


The cost function is given by equation (23) 
x2 x70) 2x, (0)x, (0) x," (0) (23) 
Tf mm ds Valli asi ui 7 ala ga a 
The-‘elements of the P matrix are those given in equation 
(22). The weighting matrix from equation (16) is 


(k-4)° (k= 4) (1-2.8) 
B= (29) 
(1-2.8)(k-4) (1-2.8)* 
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then; J~becomes 


J = x,~(0)(0.5k"-0,08k-3.248) + 2x, (0)x2(0)(0.5k-448) 4 


x, (0) (0. 5k-2. 3848) (30) 


Only points on a unit circle in the initial condition 
plane need be investigated. H#quation (23) is symmetric with 
respect to a diameter on the initial condition unit circle. 
As a further simplification, only the upper semicircle will 
be™eensa dered, 

The root locus for this system is shown in Figure 7. 
For each value of 6, where @ is the orientation angle of the 
state vector around the upper unit semicircle of the initial 
conditions plane, Figure 7(a), a value of K opt is found by 
minimizing the cost function with respect to this parameter. 
The plot of optimum regions for the corresponding poles of 
the actual system is shown in Figure 7(b) together with the 
desired model poles. The extremes are located at (-0. 5+ 
j1.936) and (-0.5+j0.775) for the upper half of the s-plane. 
Rach pair of conjugate poles on this region will give an 
optimum performance only when the system is excited by the 
initial conditions that yielded that pair of poles. For 
any other initial state there will be a different set of 
optimum poles and a correspondingly different settingser 
the free parameter for optimum response. It is not possible, 
of course, to modify the adjustment of the free parameter 
according to the type of disturbance or even worse, to pre- 


dict its values. It is then required to adopt a single 
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Fig. 6. A constrained second order system 
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Fig. 7. (a) Initial condition plane, (b) Root locus 
for s- +s @ b= 0 
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value for the free parameter that will give the best average 
index for all possible situations. Figure 8 shows the value 
of the minimum cost. as a function of initial conditions. 

Tne horizontal axis represents the orientation of the initial 
state vector. This curve shows that a vector orientation 
eur LYON w42l give an absolute minimum for the cost function 
of this example. This corresponds to Kopyt = 0.85 and the 
optimum pole lies right at the bottom of the allowable 
region. This would seem the best adjustment to use in the 
actual system. Figure 9 shows tne cost incurred by the 
system as a function of the initial state for different 
Optimum settings of tne free parameter. It is seen that 

the optimum parameter obtained at the absolute minimum of 
tne cost function of Figure 8 yields the highest index 

when used with all initial conditions. The area under 

each curve is a measure of the cost accumulated by the 
system when excited by all initial conditions on the upper 
semicircle of Figure 7(a) for a given optimum adjustment 

of the free parameter. Figure 10 shows a plot of the 

total cost as a function of the free parameter, 

The curve in Pigure 8 also exhibits another local 
minimum, which for some models is the absolute minimum of 
the cost function. The corresponding optimum parameter 
also produces a high average index as snown in Figure 9. 

The minimum area is obtained for a Kopt = 2969, Weich 


fo) 
corresponds to an initial state vector orientation of 25. 
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Fig. 38. Minimum cost versus initial conditions 
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Fig. 9. Cost as a function of the initial state for diffe- 
rent optimum settings of the free parameters 
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This particular value of k would be the design selection 
VidehMetiteae  CxiscLoeasunat Ormeprobabidaty on _thesinitial 


state 


p(0) = £ ao , lps gl 2) ap YO 


i 
1" 
For any given probability distribution on the initial 


state, p(@), we define a modified index of performance, 


which will be independent of the initial conditions, 


T 


The minimization of tnis tndex with mespect to the sires 
parameter yields an adjustment for which the average cost 
of the system is a minimum. As a consequence, the maximum 
value of the cost for any possible initial state will also 
have been minimized. This is equivalent to a minimax 
solution in optimal conexoell theoryl3, 

For the case of tne uniform prooadvility considered 


above, the modified index is 


Jae a J (6, 8) a6 


rT 
2, 
+f i cos 8 + op, Sm 8 Cos @ + Pe 2" O | de 


5 Stata (22) (32) 


The minimization of equation (32) with respect to k yields 
the result previously obtained. 
The size of the optimum region depends on the prox- 


imity of the model poles to the root locus of the actual 


Be 


system. Figure i1 shows now these regions increase as the 
model poles get further apart from the root locus of the 
actual systen. 

Only the upper half of the s-plane is shown. For 
model poles on the root locus of the actual system the 
regions are simply reduced to a point, The growth of tne 
regions as the model poles get further apart from the actual 
system root locus is more accentuated in the positive jw 
direction, an indication of the tendency of the index to 
make the natural frequency , w, » of the actual system equal 
to that of the model. The effect of the value of the 
imaginary part of the model poles on the size of the 
region is not pronounced, except that for small values 
ofw, some optimum poles could lie on the real axis. only 
near model poles are of interest nere since remote poles 
would probably require a change in the configuration of 
tne actual system, Which is outside the scope of this 
discussion, 

Tne results stated so far are all based on the minimi- 
Zatiodw of B performance index which is a function of the 
free parameter for a given set of initial conditions and 
can be represented by a two-dimensional curve. This curve 
is actually the intersection of a plane and a more general 
cost surface obtained when both parameters in the actual 
system are mee Lit If in equation (23) the actual system 


parameters are variable the cost function becomes 
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Fig. 11. Size of optimum regions for several models 
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(33) 


The shape of this cost surface will depend on the 
initial state used, as shown in equation (33) but the mini- 
mum value of this cost will always be zero and lie on the 
point (4,2.8), which corresponds to the model of this 
example. The cost surface is shown in Figures 12 and 13 
for two different values of the initial state vector orien- 
tation. It can be seen that while the absolute minimum of 
tne surface is fixed, the local minimum of the inter- 
ception of the plane a, = 1 with these surfaces changes 


2 
With initial conditions. 


THIHD ORDER SYSTEM 


Consider a tnird order linear control process charac- 


terized by the differential equation 


ef cx 2. 5X + a,x = 0 . (34) 


where the coefficient 2. is constrained by the plant param- 


eters to a value of 2. The specifications suggest a model 
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Fige 12. Isometric projection of cost surface for a model 
characteristic equation s¢ + 2.8s + 4 = 0 and an 
initial state vector x'(0) = (1520) 
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Fig. 13. Isometric projection of cost surface for a model 
characteristic equation s° + 2.88 + 4 = O and an 
initial state vector x'(0) = (0.866 0.5) 


It 


Priigeng (os 
ye 7 uy ee 4 207 = 0 (35) 


Because of the constraint in the actual system, these 
characteristics cannot be duplicated. The free parameters 
have to be adjusted to obtain the optimum performance. The 


cost function from equation (18) is 
JX “ (0): peex, (0O)x.(0) pram x, (0) 
8 = Ey WTP eee 2° 72> Me’ 3° > ye 


ee. 2 
Ka (0) pap ex, (0)x,(0)p, emo )P.. (36) 


fi elag le Pay is a particular element of the matrix P obtained 


from equation aig) thl, Tne weigntingwneeriees now 


(a, -20)* (a,-20)(a5-18) (a,-20)(2-7.8) 
B= | (a-18)(a,-20) (a-18)* (a,-18)(2-7.8)] (37) 
(2-7.8)(a,-20)  (2-7.8)(a,-18) (2-7.8)° 


lam 


The cost function is then 
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woere D= Boaz Aye 

The effects of the isAitial state on t1¢ ontinua vole 
location snowld be investigated for all voines oh the Buit 
Sonere in the state space. Because of trie Syametry of 
equetion (36) with respect to a diameter, only a hemispnere 
Meed bs considered. For two desrees of freedom, the cost 
function corresponds to a tnree-dinensional cost surface. 
The shaoe of tuls subsurface wili chanzse with initial state 
and its corresosonding miniaum point will move around on the 
B4-ap plane. The sudsurface correspoaiins to a coerficient 
az = 7.3 will have.a ainimum value of Zero at tne model | 
values, Llodedendently of tne initial state. A tyoical sub- 


surface is shown in Figure 14 for the case x4(0) = (-0.754 


0.1335 0.643). By computer searching techniques the optimnun 
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for tnis system sere found for several initial 
conditions om tie unit nenisphere to the right of tne xg 
olane. | 

The roots obtained from tnese ontinmum characteristic 
equations constlomerate in Limited regions on the s-plane. 
Tnese resioas are snom in Figure 15. Observe that the 
prooi2a of stability nas been’ autonatically solved. If the 
constraint on the actual system were released these opti- 
aum rezions Jould collapse into points where the nodel poles 


are and the cost function would be identically zero for all 
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Fig. 14. Isometric projection of cost subsurface for a mo- 
del characteristic equation s2+7.8s“+18s+20=0 and 
an actual system with a constraint a,=2. The i- 
nitial state vector is x!(0)=(-0.754 0.133 0.643) 
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Fig. 15. Optimum pole regions for third order system with 
constraint in az = 2 
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Iie lal conditions. 

Fach set of optimum poles is associated with a value 
of the cost function, wnich is a minimum for the initial 
state from which this set was obtained. For any other 
initial Seeve there will be a different set of optimum 
parameters. The final selection of parameters is based on 
tne statistical characteristics of the disturbance. The 
modified index is obtained from an extension of equation 
(31) to a three-dimensional space. For a uniform proba- 


bility distriypution over the memismnere this index is 


TW Ai 
Jmod = ae J(@, A, %) dod¢ 
a a (Pu + Pee + & Psa) + & pis (39) 


This modified index@is Independeaveor swine initial 
state, The minimum average cost function, and hence the 
Minimax solution is obtained for this example with the 


optimum coefficients 
a4 = 6.5 


aod = Se 7 


The corresponding poles are shown in Figure 15. 
A reasonable choice is to assume an initial state 


vector 
(= Glo 


The optimum set of poles for this case is shown in Figure 
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15. The cost subsurface is shown in Figure 16. This set 
of initial conditions resembles the widely used step 
function in control systems design. 

A worthwhile observation is that the available free- 
dom in the system parameters could be used to reproduce 
exactly some of the corresponding model coefficients. CGon- 
Sider the actual system and model differential equations 


as given in equations (34) and (35) 


x + 2x = aX + a,X = O (34) 


y 47. 89 2 dey +» 20; = @ (35) 


The direct reproduction of the model coefficients in the 


actual system will yield 
xX + 2x + 18x + 20x = 0 (40) 


Tne corresponding poles are close to the optimum regions 
obtained by the previous method, They are shown in Figure 
15. This possibility in parameter selection is limited by 


the stability requirements a,>0O, Ao DOs A5A3-a, >O. =when- 


1 
ever possible, this criterion gives a good, time-saving 
procedure if the thiag requitgement for stability ks weil 
mec. 

It is interesting to see how the optimum poles migrate 
as the initial state vector sweeps the normalized state 


Soace. A Mapping of a unit Cireie om the X4~X plane in 


a 
the state space onto the s-plane is shown in Figure 17 for 
this example. The optimum poles return to the same point 


after the initial conditions vector turns 180°. 
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Fig. 16. Isometric projection of cost subsurface. Model 
and actual system as in Figure 14, Initial state 
vector x4(0) = (1 0 0) 
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Consider now tnat for tne same performance specifica- 
tions, the constraint on the actual system is somewhat less 
stringent. Assume that the system's cnaracteristic equation 


is 
s? HO gES AaoS + a, = 0 (41) 


Tne model poles willl now bewe@loser to the attainarcm 
regions of the given system and consequently, the results@om 
the optimization procedure will beminecluded in smaller a= 
gions. These regions are shown in Figure 18. There is also 
a net displacement of these régiows vo tae left and nowerae 
complex comjugate model poleswapemeiemmGea, The setae 
poles obtained by tne mere duplication of model coefficiemes 
in the free terms is also shown. 

Using tne modified index of equation (39) for a uniform 


pm@bability on the initial Stave. ene wi inales-Pection gs 


a = ie 3 


a Sr oe 


Tne optimum poles are shown in Figure 18, 
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Fig. 18. Optimum pole regions for third order system 
with constraint in az = 6 
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CHAPTER IV 


WHVALUATION OF RESULTS 


Because of the constraints on the actual system and 
tne requirements for stability, the poles of the systen 
transfer function are restricted to certain regions Gmeene 
s-plane. Tne optimum pole regions nave a size, shape, and 
position according to the proximity of the model poles to 
the allowable regions. Tne selection of system parameters 
for optimum design was based on tne statistical properties 
of the disturbing initial conditions, That set of pomes 
from the optimum regions wnich gives the minimum average 
index under all possible initial conditions is the dé@simeed 
solution to the problem. 

The index of performance developed in Chapter II was 
given the form of equation (1) for mathematical convenience. 
The two indices differ in the description of the weighting 
matrix. However, the optimum parameters obtained are 
identical for every initial state. The elements of the Q 
matrix in equation (1) were chosen to place the absolute 
minimum of the cost function at the coordinates of the 
Ceetricients of The model Characterisric equationL? | Whisie 
the elements of the B matrix in equation (17) came directly 
from the general expression of the performance index of 
equation (12) when expressed in matrix form for mathematical 
manipulation, 

A comparison between the two indices is done for the 


second order system. ‘The weighting matrices, ft] and B, 


L182 


ae 
4 
1t+AQjy AQ, a2 
AQ e ot Aa, 


wheme A = 1/9? and XK, = (9° - 29) /9? 


B (94 a (94 ee) aay 
(m-22N\gu-%1) (a ~ 2) 


The expanded form of the cost function given in equation 


(2) is, by use of equatiien (22) 
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VA 2. 2, 
= x*(0) oe Be, a. 4, eke 
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- Were J. sz index of performance with matrix Q 


J index of performance with matrix B 
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Equations (44) and (45) are related by 


Ra es etc eM amon: (O amee “ (0) (46) 
See 2) 0 (7 ee > 6. /E) 


For any given initiad, conditions, the right hand sidevof 
equation (46) is a constant related to the model coeffic- 
ients and the particular set of initial conditions. Thus, 
in general, botn cost fumctionsediffer by a constant factors 
wnich is the sciitee of the model coefficient Bs and an 
additive constant given by equation (46). As a conse- 
quence, the optimum coefficients of the actual system 
obtained by the minimization of these two indices are 
identical. This is true for systems of any order. How- 
ever, the algorithm to determine the B matrix is simpler 


Charm genat fOR@ume ©) Matrix: 
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